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Abstract 

The analytic properties of a class of Generalised Husimi Functions 
are discussed, with particular reference to the problem of state re- 
construction. The class consists of the subset of Wodkiewicz's opera- 
tional probability distributions for which the filter reference state is a 
squeezed vacuum state. The fact that the function is analytic means 
that perfectly precise knowledge of its values over any small region of 
phase space provides enough information to reconstruct the density 
matrix. If, however, one only has imprecise knowledge of its values, 
then the amplification of statistical errors which occurs when one 
attempts to carry out the continuation seriously limits the amount 
of information which can be extracted. To take account of this fact 
a distinction is made between explicate, or experimentally accessi- 
ble information, and information which is only present in implicate, 
experimentally inaccessible form. It is shown that an explicate de- 
scription of various aspects of the system can be found localised on 
various 2 real dimensional surfaces in complexified phase space. In 
particular, the continuation of the function to the purely imaginary 
part of complexified phase space provides an explicate description of 
the Wigncr function. 
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1. Introduction 

The purpose of this paper is to discuss the analytic properties of a class of 
generalised Husimi functions, and the bearing that these properties have on the 
problem of state reconstruction. 

The functions we consider may be written 
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Q\e{x,p) = J dx'dp' exp 



W{x',p') (1) 



where W is the Wigner function, and where the notation xe,P9 means 

^xg\ _ f cos 9 sm9\ /x\ , , 

,) ^ y-sine cose) \p) '^^^ 

The functions Q\e belong to Wodkiewicz's class of operational probability dis- 
tributions 1^, They are the distributions which result when the filter reference 
state is chosen to be an arbitrarily squeezed vacuum state . They have previously 
been discussed by Halliwell by Wiinsche and by Wiinsche and Buzek 

In an earlier paper Q we discussed the physical interpretation of these functions. 
We showed that Q\g describes the probability distribution of measured values in the 
case when one makes a simultaneous, retrodictively optimal measurement (using 
homodyne detection for example) of xg, pe to accuracies "'""TfX ''^^P^'^" 

tively [xg, pg being defined in terms of x, p by the obvious analogue of Eq. (g)]. This 
is a universal property: it only depends on the measurement being retrodictively 
optimal, and is otherwise independent of the details of the particular measuring 
process employed. 

It is natural to ask: how much does one learn about the underlying quantum 
state from a knowledge of Qas? Or, to put it another way: what is the informa- 
tion content of Qxg7 The answer to this question depends on how accurately the 
function is known. If it is known with perfect accuracy, then it provides all the in- 
formation necessary to reconstruct the density matrix ||^, ^, 14, p^. 



However, the process of reconstruction is extremely sensitive to small variations in 
Q\g . Consequently, if one only has an approximate, experimental knowledge of 
the function, then it only provides a very limited amount of information regarding 
the state. One may express this by saying, that although the transformation by 
which one obtains Qxg from the density matrix is invertible in point of mathematical 
principle, it is effectively non-invertible in point of experimental practice. 

These statements apply to Q\g{x,p) regarded as a function of two real variables. 
However, Qxg{x,p) is in fact an analytic function of x, p which continues to a 
holomorphic function defined on the whole of complexified phase space |l^, 
[1^ , Denote the continuation Q\g{x,p). Q\g{x,p) can itself be measured 

directly (by the method of quantum state sampling, for instance [^). The question 
we address in the following is: how much information does one acquire by such 
measurements, either of the function as a whole, or of its restriction to various 
surfaces in complexified phase space? 

Of course, if one has perfectly accurate knowledge of Qxg{x,p) on any small 
region of complexified phase space having at least one accumulation point then, 
because the function is holomorphic, that suffices to determine the function every- 
where, and it consequently suffices to determine the density matrix. However, from 
our point of view this fact is not very interesting — for the reasons explained above. 
We are primarily interested, not in the information which can be extracted by some 
God-like being, endowed with perfectly accurate knowledge regarding the values of 
the function in some region, but rather with the information which is available 
to a real human experimenter, whose knowledge is necessarily imperfect. We are 
accordingly led to make a distinction between two different kinds of information. 
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which we refer to as "exphcate" and "imphcate" (this terminology is based on the 
discussion in Bohm and Hiley [|8j). Exphcate information is information which is 
experimentally accessible. Implicate information is information which, though in 
some sense present, lies hidden in the fine details of the analytic structure, and is 
therefore only accessible if one has a perfectly accurate knowledge of the function. 
In these terms the question posed in the last paragraph may be rephrased: what 
is the explicate information content of the function Q'^^g(x,p)7 And: how is this 
information localised? 

We begin in Section || by deriving what we shall refer to as the fundamental 
inversion formula. In the case 9 = this formula was obtained by Mehta and 
Sudarshan However, it appears to us that it has not received the attention 

it deserves. There do, of course, exist a number of other formulae expressing the 
density matrix in terms of the generalised Husimi function |[ Q ^ |l^, |2y, |2^, 
|2^ . However, the formula derived in Section ^ differs from all of these in that it 
involves neither an integral nor an infinite sum. The formula shows, that modulo 
a Gaussian factor, the squeezed state density matrix elements (diagonal and off- 
diagonal) may actually be identified with the function Q\g- It follows, that having 
an approximate knowledge of the values of Q1g{x,p) for x,p arbitrary complex, is 
equivalent to knowing the density matrix, up to the same degree of approximation. 
In other words, Q^g (unlike Qxe) provides a fully explicate description of the density 
matrix. 

In Section |3| we consider what happens if one attempts to use the analyticity to 
infer the complete function Q^g starting from a knowledge of the restricted function 
Qxe- We show that the errors grow exponentially fast as one moves away from the 
real plane. This means, that from an experimental point of view, the values which 
takes well away from the real plane are effectively independent of the values 
which it takes on the real plane. In Section || we go on to give a more precise 
definition of the concepts of explicate and implicate information. 

The explicate information content of Q\g is not spread out evenly (so to speak). 
Rather, an explicate description of various different aspects of the system is to be 
found concentrated on various 2 real dimensional subspaces of the 4 real dimensional 
complexified phase space. In Sections we discuss this phenomenon. In Section ^ 
we consider the problem of expressing the x^-representation density matrix elements 
in terms of Q\g- In Section ^ we consider the problem of expressing the function 
Qxe given by one choice of the parameters A, 9 in terms of the function Q\,g, 
given by some other choice. In Section ^ we consider the problem of expressing 
other phase space distributions in terms of Q\g- In particular, we show that the 
continuation to the purely imaginary part of complexified phase space contains an 
explicate description of the Wigner function. We also show that the existence and 
regularity of these distributions depends on the growth of Q^g as one moves away 
from the real plane. 

Finally, in Section ^ we show how the properties of Qxe can be used to give an 
interesting illustration of the concept of complementarity. 



2. The Fundamental Inversion Formula 
We begin by fixing some notation. Let axe be the annihilation operator defined 



by 
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'xe - yj^e + iXpt 



where xg, pg are the operators defined by the obvious analogue of Eq. (g). Let \n) 
be a normalised eigenstate of the number operator a\Qa^g with eigenvalue n, and 



3 



let \{x,p)^g) be the squeezed state defined by 

l(2;,P)Ae) =exp(-i|z|^) (3) 



where = {j^s + i^Pe)- Then the generalised Husimi function can be written 

Q\e{x,p) = ^{{x,p)xe\p\i^^P)xe) 
where p is the density matrix. More generally one may define, for any operator A, 
Axe{x,p) ^ {{x,p)xg\A\{x,p)xg) (4) 

We will refer to this function as the generalised Husimi transform of A. 

The problem of reconstructing the density matrix from an exact knowledge of 
Qxe, and more generally the problem of reconstructing an operator A from an 
exact knowledge of A\0 , comes down to the problem of using the diagonal matrix 
elements in a squeezed state representation to reconstruct the off-diagonal elements. 
The possibility of such a reconstruction may, at first sight, seem rather surprising. 
It is certainly not guaranteed just by the over-completeness of the set of squeezed 
states (one can easily think of over-complete sets which do not have this property). 
What makes it possible is the analyticity of the transform. 

There are, in fact, two different kinds of analyticity which are relevant. In the 
case of a pure state Prugovecki ||l^ has used the analyticity of the Fock-Bargmann 
wave function Q to show that the Husimi function completely determines the 
state. However, it is not clear how to extend Prugovecki's argument to deal with 
the case of an arbitrary density matrix. Also, Prugovecki does not give a formula 
for actually carrying out the inversion. 

Rather more useful is the result obtained by Mehta and Sudarshan ||l^. Since 
Mehta and Sudarshan do not consider the case 6* 7^ 0, and since we feel that the 
result deserves to be more widely known, we give the argument here. 

In view of Eqs. (||) and (^) we have 

AAe(a;,p) =exp(-|z|2) ^ ^_ Jn\ A\m) 

XJntil now we have been assuming that x and p are both real. However, we can still 
make sense of this formula for complex x, p. Define 

z± = ^ (^jX0 ± iXpe^ (5) 

If X, p are allowed to be complex then z± are independent complex variables. 
Consider the expression 

^Ae(a;,p) =exp(-z_z+) V {n\A\m) (6) 

— ' y/nl m! 

n,m— 

It is straightforward to show that the series on the right hand side of this 
equation converges if there exists a positive constant K, and a constant a in the 
range < a < 1, such that 

\\A\n)J\<{Kn'^r 

for all n. In that case Eq. ^ gives the continuation of A\g to a holomorphic 
function defined on the whole of C^. 
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Let us now see how we can use the continuation to invert the transform. Consider 
the off-diagonal matrix element ((a;i,pi)^g| A \{x2,P2)xe)- We have 

{ixi,Pi)xe\A\{x2,P2)x0) 



' vn'.m] 

n,m— U 



where 



^ (^jXr0 + iXpre^ r = l,2 

The points (a;i,pi) and {x2,P2) both belong to M^. We now use them to define a 
single point {x,p) in complexified phase space 

X — —{xi + X2) + — sinh 2ri sin 29 {xi — X2) — — (cosh 2ri + sinh 2ry cos 26) (pi — p2 ) 

1 z ^ 

p = 2(f^i +-P2) + 2 (cosh27y — sinh 2?7 cos 26') [xi — X2) — - sinh 277 sin 26* {pi — P2) 

where we have set 77 = In A. We have 

1 iX^ 
xe = COS0X + sin6'p = - (xw + X20) ^ {pie - P2e) 



pg = - sine x + cos0p = i {pig +p2e) + ^ {xw - X2e) 



Consequently 

z_i_ =^ Z2 and z_ = z* 

where z± are the quantities defined by Eq. (^). Inserting these relations in Eq. 
and comparing with Eq. (^) we deduce 

{ixi,Pi)x0\A\{x2,P2)xe) = exp [-i (|zip + |z2p) +21^2] ^Ae(a^>P) 
which can alternatively be written 



{ixi,Pl)x9\ 


\A\{x2,P2)xe) 




\{x2,P2)xe) 



A\e{x,P) = —r, tSt 



Eq. (H) relates the off-diagonal matrix elements to the analytic continuation of the 
diagonal elements. We will refer to it as the fundamental inversion formula. 

If 61 = and A = 1, and if A is a density matrix, then the expression on the left 
hand side is the Q-function, continued to complex values of x and p (and modulo a 
factor 2^); while the expression on the right hand side is the i?- representation ]25[| 
[modulo a factor exp(— zjf Zj)]. We see, therefore, that the i?-representation is es- 
sentially the same thing as the analytic continuation of the Q-function — as was 
originally noted by Mehta and Sudarshan [|lO| . 

There do, of course, exist other ways of inverting the transform ||, ||, |l^, 
[T9I , pO| , |2l| , p2| . However, these methods all involve, either an integral, or else the 
calculation of an infinite sum. The equation just derived is more straightforward. 
It shows, that in order to invert the generalised Husimi transform, all that one has 
to do is to perform the analytic continuation, and then to multiply by a Gaussian 
factor. 

Eq. (||) also shows that the function can be measured directly, by the method 
of quantum state sampling 
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3. The Physical Significance of the Analyticity 



Eq. (g) shows that an exact specification of Qxe{x,p) for real x, p suffices to de- 
termine the density matrix. Actually, a very much stronger statement is true: once 
is exactly known on a region of complexified phase space having at least one 
accumulation point, then that suffices to determine the function everywhere, and 
it consequently suffices to determine the density matrix. Alternatively, an exact 
knowledge of the function together with all its derivatives at any one point suffices 
to completely determine it at every other point (as has been stressed by Wiinsche 
and Buzek Q). The function might therefore be compared with a hologram: infor- 
mation sufficient to reconstruct the whole picture is folded into every little piece of 
it. 

Qxe is a very intricately structured object. In particular, the analyticity means 
that there exist subtle connections between the probability of finding the system in 
one part of phase space, and the probability of finding it another. As an illustration 
of this point, specialise to the case 9 = 0, and consider the following expression O] 



for the quantity ;i^Q(n| p \n 
of the last section) : 



/AO 



'AO 



being the number state defined at the beginning 



AO 



M p m 



AO 



2tt 



^ in- r)\ (r!)2 T 

r=0 ^ I \ ) 



Qxq{x,p) 



(9) 



x—p—Q 



xo{n\ p |"-)ao probability of finding the system in the state \n)^Q. It therefore 

tells us something about the probability of finding the system in the vicinity of the 
surface ^a;^ -I- X^p^ = 2n + 1. Yet in order to calculate the probability we have 
to evaluate Qxo ^md its derivatives, not in the vicinity of this surface, but at the 
origin. 

The feature just evinced could be regarded as a kind of non-locality. It is some- 
what reminiscent of the kind of non-locality that is manifested by a violation of the 
Bell inequalities p6| . 

Bohm and Hiley ||l8| place much emphasis on what, following Bohr ||2^, they 
describe as a property of "wholeness" possessed by quantum mechanical systems. 
In their discussion of this feature Bohm and Hiley also make use of a hologram 
analogy. They argue that the property of wholeness does not essentially depend on 
assumptions peculiar to "hidden-variables" interpretations of quantum mechanics, 
such as the Bohm interpretation. On the contrary, they maintain that it appears, in 
one guise or another, whatever the interpretational scheme one adopts (it appears 
in the Copenhagen interpretation favoured by Bohr, for example). The analyticity 
of the function Q^g could be regarded as another illustration of these ideas. 

This phenomenon of "wholeness" is certainly intriguing. However, it is also 
elusive. The non-local connections in the Bohm interpretation cannot be used for 
faster-than-light signalling. The practical significance of the kind of non-locality 
exemplified by Eq. (^ is also very limited. It critically depends on Q^g being known 
with perfect exactness. However, Q^g could in fact only ever be measured up to a 
certain non-zero error, which seriously restricts the usefulness of the analyticity. 

To see this, consider the problem of calculating at a point (a; + ?7,p + C) 
using the expansion 



Ql9{x + V,P + = 



oo 

y 

n\va\ dx^dp^^ 

n.m—O 



Qxe{x,p) 



Suppose that the derivatives appearing in this sum have been measured indepen- 
dently. For the sake of example, suppose that the error in the measurement of 
Q^nQj^m Q\0 is ±A'"^"cr, where cr is a dimensionless constant. Let AQ be the error 
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in Qxsix + ?7,p + C), as calculated using the above formula. Then 



\n,m=0 ^ ' J 




where /q is a Bessel function of imaginary argument p8|] . Asymptotically ||2^ 



^Q^\o[A^\\C\) 'exp(A-iH+A|C|) 



It can be seen that the error grows exponentially as one moves away from the 
point (a;,p). If, instead of direct measurements of the derivatives, one only had 
imprecise information regarding the values of the function itself in some small region 
surrounding the point (a;,p), from which the values of the derivatives had to be 
estimated, then one would expect the error to grow even more rapidly. 

We see that the situation is not really so different from the case of a function 
which is not analytic. Given an experimental determination of the function in some 
region, one can use the analyticity to infer something about its behaviour a small 
distance outside that region. However, one cannot go very far before the data is 
swamped by the errors. 

Although they are not independent in point of mathematical principle, it may 
be said that the values which takes in well-separated regions of complexified 
phase space are effectively independent in point of experimental practice. 

The situation here could be compared with the case of a classically chaotic 
deterministic system. In principle the state of such a system at any one time 
suffices to completely determine its state at every other time. However, the practical 
usefulness of this fact is very limited. The information is in some sense present, but 
mostly not in a form which it is possible to access. 



In view of the discussion in the last section we are led to distinguish between 
two different kinds of information which, following Bohm and Hiley [T^ , we will 
refer to as "explicate" and "implicate" (although we do not use these words in quite 
the same sense as Bohm and Hiley). Explicate information is information which 
approximately survives the process of making inexact measurements. Implicate 
information is information which is experimentally inaccessible. 

In order to make these notions a little more precise, consider a physical quantity 
which is given by an expression of the form 



In this expression the domain of integration D might be a 4 real dimensional volume 
in complexified phase space, or it might be a lower dimensional surface (such as 
the real plane). / might be an ordinary function, or it might only be defined in a 
distributional sense. 

Suppose that / is an ordinary function. In that case the error in the integral is 
of a comparable order to the error in Q\q- The same is true for certain, not too 
singular distributions. It is true if / is a ^ function, for example. In such a case 
we will say that Q\g^ restricted to the region D, contains explicate information 
regarding the physical quantity in question. 

Suppose, on the other hand, that / is a very singular distribution. Suppose, for 
example, that / involves many high order derivatives of the b function. In that case, 
although the integral may still be of some theoretical interest, it cannot easily be 
used to make practical calculations based on inexact experimental data concerning 



4. Explicate v. Implicate Information Content 
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the function Q'xg. We will say that Qxe^ restricted to the region D, only contains 
implicate information about the physical quantity in question. 

The distinction between "explicate" and "implicate" is qualitative, and it should 
not be taken in too absolute a sense. For instance, one cannot entirely rule out the 
possibility of measuring all the quantities on the right hand side of Eq. (^), and 
using them to make a good estimate of ;^o("l P I"')ao' However, it would obviously 
be very difficult if n is large. 

As an illustration of these concepts, consider the problem of calculating the 
expectation value of an operator A. Using Eq. we have 

Tr(i/5) = ^ J dxidpidx2dp2 {{xi,pi)^g\ A\{x2,P2)xe) {i^2,P2)xe\ p\i^i^Pi)\e) 
^ d^xd^pexp(^^x^,e~2\^p^,s]Als{x*,p*)Qlg{x,p) (10) 



TT 



In this expression xi, pi denote the imaginary parts of the complex variables x, 
p. xie, Pie are defined in terms of xi, pi as in Eq. (^). The integral is taken over 
the whole of complexified phase space. It can be seen that Q^g contains explicate 
information about {A) whenever is defined as an ordinary function — which is 
the case for most operators of physical interest. 

Suppose, on the other hand, that one wants to calculate {A) using only the 
values of the restricted function Qxq. Inspection of Eq. (|l^) shows, that due to the 
presence of the Gaussian factor, the integrand will be strongly peaked about the 
real plane, provided that ^^g(x*,p*) does not grow too rapidly with increasing Xj, 
Pj. In that case we may expand ^^g(a;*,p*), and Qxg{x,p) about the point {x^,p^) 
(where x^, p^ denote the real parts of x, p): 

/j2i2_ / ^2 n\22\v^/ ' 



Tr(ip) = ^ / '^'^ exp (-^xie - 2AVe 

n,n\m,,ni'—0 



n\ n'\ ml m'l 



I I fin+m Qn'+m' 

X 2;"+" , ^ Axeixr^^pr.)—, jQxeix^.PT,) 



One finds 



Tt{Ap) ^ J dxB.dpB,Alg{x^,,pB.)Qxe{xK,pn) (11) 
where A^g is the generalised anti-Husimi transform ||2^, |2^, defined by 



^xei^niPR) = exp 



_1/ 2 52 ^ I d 



2 V" dxlg dpi 



Ax9{xT^,Pn) (12) 



There are, of course, easier ways [|| of deriving Eq. (|ll|). The advantage of the 
derivation just indicated is that it brings out the connection between Eqs. ( p^ ) 
and (pT|). It also brings out the fact that the existence of A^g as an ordinary 
function depends on the growth of as one moves away from the real plane. We 
will return to this point in Section 0. 

If A'^g exists as an ordinary function, or as a not too singular distribution, then 
Qxe provides us with explicate information regarding (A). However, A^g is in fact 
often very singular. The conclusion is, therefore, that whereas Qas provides us with 
an explicate description of the quantum state, the description provided by Qxe is, 
for the most part, only implicate. 

Qxe is the probability distribution describing the outcome of a retrodictively 
optimal determination of phase space location 0| . It might therefore be compared 
with the function | {x\ Tp) p. The values which takes away from the real plane 
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provide us with the additional, exphcate information needed to calculate the ex- 
pectation value of an arbitrary observable. They might therefore be compared with 
the phase of the function (a; | V') • 



5. Connection with the ^^-Representation Matrix Elements 



The advantage of Eq. (|11|) is that it only requires us to integrate over two real 
variables, unlike Eq. ( [l0| ) which requires us to integrate over four. The disadvantage 
is that A'^g is often extremely singular. It is natural to ask whether one might do 
better by integrating over some other surface, instead of the real plane. In the 
next three sections we investigate this question. We begin by identifying a surface, 
having two real dimensions, which contains an explicate description of the matrix 
elements 

(where la;)^ is an eigenket of X(j, with eigenvalue x). The result we derive is a 
generalisation of the one proved by Davidovic and Lalovic jl^ . 
It is readily shown ^ 

J^x'\ix,p)^g) = {X'\{x^,p^)^g) 



^/tt (Acos(5 — zA ^sin5) 
X exp 



2 (A cos J — iA~i sin (5) 
where 6 = 6 — (p. Consequently 



A ^ cos(5 — iAsin(5 , , ,1 

[X - x^j + ip^x - -p<px^ 



Qxe{x,p) 



J dx'dy' exp -i {{x' - x^f + y'^) 



~y'{x' - x^) + 2ip^y' 



where 

b^\^ cos^ 5 + sin^ 5 
c = - (A^ - A"^) sin 5 cos ^ 

Continuing to complex values of x and p we find 

exp(— feu^)(35,6i sin + md, v cos cj) + iuf) 



l^x' -y'\p\x' + y')^ 
(13) 



7r2 \/b 



dxdy exp [2i{ux + vy)] 

X exp [-\{x'^ + 2icxy + y^)] {x -y\p\x + y) 



where 



d = h cos ( 
f — b sin (J 



' c sin (j) ~ cos cos 6 + \ ^ sin 9 sin 6 
c cos 6 — sin 6 cos 5 — A~^ cos 9 sin S 



(14) 
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Inverting the Fourier transform gives 

exp [-l{x'^ + 2icxy + y'^)] J^x -y\p\x + y) ^ 

' f 

/ dudv exp [—2i{ux + vy)] 

J J 

X exp(— 6w^)(5^g [—V sin (f> + iud,v cos (f> + iuf) (15) 

which is the desired result. 

By essentially the same argument one can derive an analogous formula relating 
the a;0-representation matrix elements of an arbitrary operator A to the function 

6. Qye' IN TERMS OF 

Using similar means one may obtain an explicate description of the function Q\e 
by integrating Q'Xig, over a surface having 2 real dimensions. However, the general 
formula is rather complicated. We will therefore confine ourselves to showing how 
Qxe may be expressed in terms of the analytic continuation of the Q-function, 
= Qio- 

Suppose that A < 1 . Eq. (|lj) gives 
exp{—X~^u^)Q'^^g{—v sin 6* + iu cos 9, v cos 9 + iu sm9) 

dxdy exp[2i(A^^ux + vy)\ exp[— A~^(x^ + ?/^)] J^x — y\p\x + y) g 



7r2 A 



while Eq. (|T|) gives 



exp 



{x^ + y^)\ lx-y\p\x^y)^ 



du'dv' exp [—2i{u'x + v'y)] exp(— w^) 



X Q'^{~v' sin 9 + iu' cos 9, v' cos 9 + iu'sin( 



Combining these formulae we obtain 



exp 



1 



1-A2 
A 

^(1-A2) 



ix,p) 



du'dv' exp 



2i 



1-A- 



:U xg 



-{v'^pgf 



X exp 



1 



1-A2 



Q'^{—v' sm9 + iu' cos 9, v' cos 6* + iu' sint 



(16) 



If A > 1 we have 
1 

exp 



1- A-2'" 

A-i 
7r(l- A-2) 



\e{x,p) 

du'dv' exp 
1 



2i 



1 



X exp 



1- A-2 

,2 



upg 



-W-Xgf 



1 - A- 



A2 - r 

tW Q'^{v' cos 9 + iu' sm9, v' sin 6* — iu' cos 9) 

(17) 

Q describes the result j7| of measuring xg , pg both to a retrodictive accuracy of 
±-^. Eq. (|l6|) shows that continuing to one part of complexified phase space gives 
a more accurate explicate description of xg. Eq. ( p^ shows that continuing to a 
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different part of complexified phase space gives a more accurate explicate descrip- 
tion of pg. An explicate description of the result of every possible measurement is 
to be found somewhere in complexified phase space. 

Analogous formulae hold for the generalised Husimi transform of an arbitrary 
operator A. 

7. The Generalised s-Transform 
The generalised s-transform of an operator A is defined by 



1 

2^ 



dx' dp' exp 



,(p^'_V) + -(AV^ + A-V^) 



TT[D,,p,A 



where Dxp — exp[i(pi — xp)] is the displacement operator. 
If s > — 1 it is straightforward to show 



Axe(x,p) 



1 



dx' dp' exp 



7r(l + s) 

This formula can be inverted using 



1 



(1 



{\-\x',~xef + \\p'e~po?) 



x4^j(x',p') 



Axe{x,p) = exp 



1 92 



4 \ dxj ^ A2 dpj 



Axeix,p) 



(18) 



(19) 



It is to be observed, however, that the expression on the right hand side of Eq. ([T^ ) 
contains derivatives of all orders, which means it is exploiting information that is 
only present in implicate form. If one wants a formula which only depends on 
information present in explicate form it is necessary to make use of the analytic 
continuation. 

Continuing to imaginary values of x and p in Eq. ( p^ ) gives, after rearranging, 

1 / 1 _V V ^ J J \ AO i ■ ■ \ 

Axe[ix-e,ip~ej 



exp 



1 + s 
1 



7r(l + s) 



{X-^x^ + aV) 
dx' dp' exp 



2i 



X exp 



+ 

1 



{\-'^xx' + X^pp') 



{X-^x'^ + AV) 



1 + s 

Inverting the Fourier transform and changing the variables of integration gives 
1 



exp 



l + s 
1 



7r(l + s) 



{X-^xl + X^pI) 
dx' dp' exp 



A^xl{x,p) 
2i 



X exp 



(1 + ^) 
1 



(A "^xex'g + X^pep'g) 



1 



{X-^x'l + X'p'D 



Alg{ix',ip') (20) 



If s = and A is the density matrix this becomes 
exp[- {X-^xI + X^pI)]W{x,p) 

= — y dx dp exp \—2i{X^'^X0x'g + }?pgp'g)\ 

X exp [- (A-242 ^ y2^,2^^^ Q\g{zx', ip') (21) 

We see that the continuation of Q\e to the purely imaginary part of complexified 
phase space contains an explicate description of the Wigner function. 
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Eq. ( pOD also implies a convenient criterion for the existence and regularity prop- 
erties of the generalised s-transform, in terms of the growth of A'^^g as one moves 



away from the real plane. Specifically: A^^g {x, p) exists as a tempered distribution if 
and only if exp — Y^(A~^Xg + X^Pg) A'^^g{ix, ip) exists as a tempered distribution 
{x and p both real). 

In view of Eq. (|^) the condition may be rephrased: the necessary and sufficient 
condition for A^^g{x,p) to exist as a tempered distribution is that 



exp 



1 



-(A 



-2 2 I \2 2 
Xo+X Vg 



{i-X,-p)x6\A\{x,p)^g) 



exists as a tempered distribution. In particular, the generalised s-transform exists 
as a tempered distribution whenever the operator A is local, in the sense that it 
takes localised wave packets into wave packets which are again localised, in the 
same region of phase space. 

8. Complementarity 

We conclude by showing how these considerations can be used to illustrate the 
concept of complementarity. 

Suppose one took a classical system characterised by the classical distribution 
r and measured the observables xg, pe to accuracies i"^, vfX '^^spectively. The 
probability distribution for the values of x, p as calculated from the measured values 
of xe , Pe would be 



pxe{x,p) = ^ J dx'dp' exp [-A" 



\x'e - xgf - X'ip'e - pef] T{x' ,p') (22) 



Comparing the defining equation for Qxe, Eq. (Q), with this equation it can be seen 
that the only difference consists in the replacement of the classical distribution T 
with the Wigner function W . It is natural to ask: can one infer that the underlying 
distribution is W rather than F from experimental determinations of the functions 



Q 



xe ■ 



The answer to this question is that the non-classical nature of the Q 



xe 



does not become experimentally manifest for any single pair of values of A and 9. 
However, it does show up if one makes measurements for many different values of 
these parameters. 

To see this, consider Eq. ( p^ ) for the case 4> = 9: 

1 



Qxe{x,p) 



— 3 — / dx'dp' exp 
7r2 A J 



A2 



{{x' - xe)'^ + y''^) + 2ipey' 



y \ P \x +y\ 



In the limit as A 



and 



0FA 



exp 



6{x' - Xe) 



^X 



exp 



A2 



y' + 2ipey' 



exp [-XVe] 5{y') 



Consequently 



Qxe{x,p) ~ ^exp [-X^pI] 



ei^e] P \xg)g 



(23) 



in the small A limit. An experimental determination of the functions g{xe\ p \xe)e 
provides one with enough information for an approximate reconstruction of the 
quantum state ^, |2^, |3^, It follows, that an experimental determination of 
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the Qxe(x,p) (with x, p both real) also provides one with enough information for 
such a reconstruction, provided that the functions are known for sufficiently small 
A, and sufficiently many different values of 6. 

As an example, consider the Fock state \n) = ki)io- The Q-function is 

Q{x,p) = Qio{x,p) = T^.E^e-"" (24) 
27rn! 

where E = ^(x^ +P^) is the classical energy. Using Eqs. ( |l6|) and (H^) one finds 

^ , , sechri Itanhnl" /, , „ ,i 
Qxeix,p) = ' , il„ |cosech277|- z^e 

X exp 



2 



-isech77(A ^xl + Xpl 



(25) 



where rj — \nX, z\g ~ "^("^ iXpo) and is a Hermite polynomial. Qio is 

peaked about the classical orbit at E — n. The behaviour of Q\g for other values 
of A, 9 is illustrated in Figures ^ and 

The functions Q\e are obtained from the Wigner function by smoothing it in 
various ways. Different smoothings extract different features of the information 
present in the Wigner function. No single function extracts all the information (in 
explicate form). If one wants a complete description, then it is necessary to consider 
the complete set of functions (or else to use the continuation to complex values of 
X and p) . 

The Wigner function describes the state much more efficiently, using the single 
function W, instead of the two-parameter family of functions Qxe (restricted to real 
values of their arguments). This is an important advantage of the Wigner function. 
However, it can also be a disadvantage since, precisely because it represents the 
information in a very compact form, the Wigner function can be hard to interpret. 
The function Q\g, by contrast, describes the result of making a retrodictively op- 
timal determination of phase space location. As such it has an immediate physical 
interpretation, and the picture it presents is much easier to assimilate intuitively. 

The functions Q\g can be thought of as a set of alternative views, or perspectives 
of a single underlying object. They might be compared with a set of photographs 
of a building taken from various directions. The description provided by any single 
photograph is incomplete. In order to get an adequate impression of the building 
as a whole, one needs to see how it looks from every direction. 

It appears to us that this discussion provides some additional insight into the 
intuition underlying Bohr's concept of complementarity. In his book Interpreting 
the Quantum World Bub j32j uses Escher's well-known print Waterfall as "a visual 
metaphor for a quantum world, in which the properties of physical systems 'fit 
together' in different classical or Boolean perspectives that cannot be put together 
into a single Boolean framework". The above considerations may be taken as an 
alternative illustration of what is, we believe, essentially the same point. Measure- 
ments performed for one particular pair of values oi X, 9 are not sufficient to exclude 
the possibility that the observed probability distribution results from an underlying 
classical distribution F via Eq. (p^. It is only when one considers the complete 
set of functions Q\e that one realises that they cannot all be accounted for in the 
same way, in terms of a single classical distribution F (except in those cases where 
the Wigner function is itself a possible classical distribution p3{). 
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Figure 1. The function Qxo for the state \n)-^Q, with n = 6 and 
9 = 0. In the top diagram A = 0.5. In the bottom diagram 
A = 0.25. See Eq. It can be seen that at A = 0.25 the 

function is beginning to take its asymptotic form, as given by 
Eq. dH). 
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Figure 2. The function Qxo for the state |?t.)]^q, with n = 6 and 
A = 0.25. In the top diagram 9 = ^. In the bottom diagram 
9 = ^. See Eq. (p5|). It can be seen that the effect of varying 9 
is simply to rotate the distribution. This is a consequence of the 
rotational symmetry of the state |n)]^o. 



